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ACOUSTIC RAYTRACING TO FIXED GROUND 
RECEIVERS (PROGRAM AIRRAY) 


INTRODUCTION 


AIRRAY is a FORTRAN-IV program which adapts some of the techniques used for some 
time im underwater acoustic ray computations to atmospheric ray tracing. The basic ray path 
integration algorithm of R. J. Thompson,':? which has been incorporated into this program, 
computes the paths of sound rays reaching the ground from a source at an arbitrary altitude, 
and traveling in a fixed vertical plane. The effect of the horizontal component of the wind in 
the source-receiver plane, at each altitude, is considered in the program; other wind com- 
ponents are neglected. The present program is designed to find all eigenrays (rays connecting 
source and receiver) reaching a set of equally-spaced receivers along the ground track, either 
directly or by multiple surface reflection, and compute their travel times, ray intensity levels, 
ray angles and number of surface reflections. Eigenrays are found, and their trajectory parame- 
ters are computed, by linear interpolation in the set of rays traced to the ground. After the 
eigenrays along a user-specified set of tracks have been computed, a table of CW transmission 
losses vs. range and bearing is computed by random-phase addition; or, for a given bearing, a 


table of transmission loss vs. range and source altitude may be computed. 


Since this program will often be used to model explosive or moving sources of sound, the 
eigenray listing may be a more useful output than the transmission loss table. The transmission 


loss calculation is strictly valid only for a stationary continuous source, however, the table may 
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still be useful in indicating the locations of shadow zones on the ground. Transmission losses 


greater than 120 dB are not printed, so that high-loss regions stand out in the table as blank 
areas. 
The physical input data required is a weather profile, i.c., a table of temperature, wind 


speed, and wind direction as functions of altitude. The receivers are always on the ground, and 


the ground is flat, although it may be at an arbitrary elevation with respect to sea level, where 


the weather profile is assumed to start. 


THEORY 


An outline of the theory will be found in Refs. 1 and 2. It is of interest here, however, to 
show a comparison between the equations used here and those of single-profile underwater 


acoustic ray tracing with currents neglected, where Snell’s law applies. The basic differential 


equations in the Snell’s law case are :? 


dx 

oe (1A) 
oer ae 5 pees (1B) 
dt c?2 


and an invariant quantity along the trajectory of a ray is 


c 
cos 6 (2) 


Cm 


In the present model of an atmosphere with wind, the differential equations are 


dx — u(l + ku) — ke? 

dt 1+ ku (3A) 
dz eV (0 + ku)? — kc? 
Wen 1+ ku Ge) 


ee ce aa ee ae ae eT RE RON Te Nee ee OS, 


NRL MEMORANDUM REPORT 4185 


and the invariant quantity along a ray trajectory is 


2 


1 c 
inn (4) 
cy? Ve? + (c? = w)tan? 6 


In all of the equations above, x and z are the horizontal and vertical coordinates, respectively, 
c = c(z) is the soundspeed, uv = u(z) is the in-plane component of wind speed, and @ is the 
ray angle with the horizontal. Note that the present model is not the small-angle model some- 
times used, where Eqs. 1 and 2 are used with the still-air sound speed c being replaced by c + u 
at every altitude in the vertical plane.4 The present model treats the wind component in the 


plane exactly (Eqs. 3A and 3B), although neglecting the ordinarily smaller effect of crosswind. 


The interpolating function used here is piecewise-linear; i.e., the ray path equations are 
solved in closed form in each profile layer, under the assumption that c and u vary linearly with 


z within the layer, and are continuous at layer boundaries. 


INPUTS AND OUTPUTS 
The input card deck structure for AIRRAY is as follows: 


Card 1: Any alphanumeric title in cols 1-72. 


Card 2: DR, R2, (H(I), 1 = 1, 8) FORMAT (10F8) 
DR is the initial receiver range (and range increment) in km. 
R2 is the final range in km. 
H is the array of source altitudes in m. 
Additional source altitudes may be entered following card 2 in format (10F8), up to 
a total of 30. Always terminate the array of source altitudes with a blank or nega- 
tive number, using an extra card if necessary. Note that the source altitudes 


entered are measured from ground level. 
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Card 3: 


Card 3A: 
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(B (I), I = 1, 16) FORMAT (16F5) 


(B (ID, I = 17, 30) FORMAT (16F5) 
B is the array of track bearings (source to receivers). Use card 3A if more that 16 
are to be computed. Always terminate the bearing input with a negative number 


unless the full 30 are entered. 


(Y ®,T (@, DW, W @,1 = 1, 50) FORMAT (4F8) 


Y is the array of weather profile altitudes in m. (Y is measured from sea level—see 
below.) 
T is the array of temperatures in degrees C. 


D is the array of wind directions (bearings from which the wind blows) in degrees. 


W is the array of wind speeds in knots. 


The weather profile is read in, one altitude per card, from sea level up. The 
last altitude should be above 100,000 m, and its temperature should be the max- 
imum on the profile. Always terminate the profile input with a blank card, unless 
all 50 altitudes are used. Since the first weather profile card is for sea level, the 
first profile altitude is always 0. If a non-zero ground elevation is desired, enter it 
in m. in the Y-position on the first profile card. The rest of the altitudes in the 
profile will be reduced to altitudes above ground level in the program. Source alti- 


tudes are always measured from ground level when entered. 


The program finds and lists all eigenrays for a track (source-receiver bearing) before going 


on to the next track. (The listing for each eigenray includes the number of ground reflections 


the ray has had up to the receiver point, and the ray’s source and receiver angles, travel time, 


intensity level, and mean horizontal speed.) After all tracks for a given source altitude are 
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done, a table of transmission loss in dB vs. range and bearing is computed by random-phase 
addition of eigenrays at each receiver, and printed. Then the program proceeds to the next 


source altitude and repeats the computations for the same set of bearings. 


If only one bearing is entered, transmission loss is tabulated vs. range and source altitude 


rather than range and bearing. In this case the table is not printed until the very end. 


Appendix A is a listing of the program; Appendix B is a listing of a sample data deck with 


the resulting output. 
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